NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



34902 -H002-R0-00 


A I 


1 THDRi - 


I 



ADVANCED FLIGHT DESIGN SYSTEMS 
SUBSYSTEM PERFORMANCE MODELS 


JUNE 1980 CONTRACT NO. NAS9-15793 

(N ASA-CR-1 607b 9) ADVANCED FLIGHT DESIGN N80-31045 

SYSTEMS SUBSYSTEM PERFORMANCE MODELS. USER 
GUIDE: ENVIRONMENTAL ANALYSIS ROUTINE 

LIBRARY (IRa Defease and Space Systeas Unclas 

Group) 138 p HC A07/MF A01 CSCL 0 oK G3/54 28522 

USER GUIDE 

ENVIRONMENTAL ANALYSIS ROUTINE LIBRARY 


Prepared By 
K. C Parker 
J. O. Torian 

Systems Engineering and Analysis 
Department 


o***mu amd smc* rrsrmm amxm 




34902 -H002-R0-00 




ADVANCED FLIGHT DESIGN SYSTEMS 
SUBSYSTEM PERFORMANCE MODELS 


JUNE 1980 CONTRACT NO. NAS9-15793 


USER GUIDE 

ENVIRONMENTAL ANALYSIS ROUTINE LIBRARY 


Prepared By 
K. C. Parker 
J. O. Torian 

Systems Engineering and Analysis 
Department 

TRW 

OfffNSi AND SMCI SYSTfMS G*OUP 




L 


PREFACE 


Subsystem performance analysis Is required in Fllqht Design to assess 
the capability of the Environmental Control and Life Support System (ECLSS) 
to support the flight requirements and define operational procedures under 
contingency flight conditions. Current ECLSS modeling techniques are limited 
in the variety of configurations and they employ batch mode computer programs 
execution methods. Future spacecraft will require analysis of both a greater 
variety and a greater number of ECLSS than fur previous spacecraft programs. 
Improvements in the variety of configurations that can be modeled and a reduc^ 
tion in effort required for modeling and analysis can be accomplished by 
developing a modular computer program which operates interactively. 

An effort has been conducted to develop a modular interactive ECLSS per- 
formance analysis tool. The final reports on the effort are included in an 
Executive Summary and two Technical Reports. The Technical Reports include a 
User Guide and a sample model. 

The Executive Sumnary presents an overview of the effort. 

This Technical Reports presents a User Guide which, due to the modular 
nature of the Program Library, includes a greater degree of technical detail 
than one for a conventional program. A sample model report supplements the 
User Guide and illustrates a complete ECLSS model set up and execution. 
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1. INTRODUCTION 




This report presents a user guide for a library of Interactive computer 
routines used to develop performance analysis models cf specific Environmental 
Control and Life Support Subsystems (ECLSS). This volume Is supplemented by 
a second volume in this series of reports which presents an example of a 
complete model set up and execution. 

The Environmental Analysis Routine Library (EARL) is designed such that 
additional ECLSS component performance routines may be added as required. To 
facilitate report revisions associtated with such additions a page. Table, and 
Figure numbering system based on the sections numbers is used. In this system 
the last number identifies the page. Table, or Figure number for that section 
whose number preceeds. This system facilitates revisions and additions with- 
out requiring complete renumbering. 


2. FORMULATION OF ENVIRONMENTAL CONTROL AND 
LIFE SUPPORT SYSTEM ANALYTICAL MODELS 


Evaluation of various Environmental Control and Life Support Systems {ECLSS ) 
performance may be conducted by the application of the subject Interactive com- 
puter program with which the user accesses a library of routines simulating the 
performance of various components and functions common to ECLSS. These routines 
are assembled with a driver routine (MAIN) to simulate the particular ECLSS under 
consideration. The assembled program Is then loaded and executed to produce 
the transient performance parameters of the ECLSS under prescribed boundary 
conditions. The contents of the MAIN are typified on Table 2.1. 

The assembly procedures for a program are shown in Figure 2.1. The master 
library of routines is extracted from a secure file. The user has the option 
to enter a MAIN routine (as for initial development of an ECLSS model), or 
extract a particular MAIN from his individual library (as for update/edit and/ 
or additional studies with a previously developed ECLSS model). The extracted 
MAIN may be altered as part of the update/edit process. The program is then 
MAP'ed and the MAIN may be stored in the user's file for future use. The 
particular ECLSS program is then ready for execution. 

The execution procedure including a variety of input/output options is 
shown in Figure 2.2. The component characteristic data and initial conditions 
may be read in from restart data stored previously or entered directly. If the 
user desires, the system will output a schematic of the ECLSS modeled. The 
user has the option to select particular nedes (component locations) to be 
included in tabular output or the system will default to include output for all 
nodes defined i the model. If plots are desired, the user simply defines the 
particular parameters to be plotted. Restart data may be stored for future 
use. Up to this point the program is executed in an interactive mode. The 


program then transfers to a second stage of execution. 

The second stage of execution Is passive In the sense the data Is processed 
with no Interaction on the part of the user except to produce hard copy of the 
output. The processing accesses Electrical Power System (EPS) data and/or 
trajectory data automatically, If required. The data can be accessed from tape, 
secure files, or Interface to other programs resident in the system. This stage 
of execution produces the tabular and plot data. 


Table 2.1. Typical MAIN Contents 


CALL START 
333 CONTINUE 

CALL STEP( . . . . ) 

3133 CONTINUE 

CALL LOOP( . . . . ) 

CALL PLATE( . . . . ) 

t « « * 

CALL MIX( . . . . ) 

CALL CONVRG( . . . . ) 

CALL PRINT 
GO TO 333 
END 


} PROGRAM CONTROL 

} CALL TO INPUT UTILITY ROUTINES 

j CONVERGENCE CONTROL 

CALL TO COMPONENT ROUTINES 

} CONVERGENCE CONTROL 
} PROGRAM CONTROL AND TIMING UPDATE 


2.3 
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Figure 2.2. ECLSS Program Execution Procedure 
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3. LIBRARY OP ROUTINES 


TWo catagorles of Routines make up the library. The first category Is 
Referenced Routines, which are directly referenced In the user's WIN. The 
second category Includes the Unreferenced Routines which are automatically 
executed, but are not referenced directly by the user. All Referenced Routines 
are discussed In the following section. Those Unreferenced Routines with which 
the user may desire Indirect communl cation are discussed In the subsequent 
section. Unreferenced Routines also Include a variety of Control, Computational, 
and Support Routines. 
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3.1 


3.1 REFERENCED ROUTINES 


This section presents a description of the functions, application, and 
parameters associated with those Library Routines which are directly referenced 
In the user's MAIN. These Include Program Control Routines, Component Perfor- 
mance Routines, and Input Utility Routines. The following discussion Includes 
the reference procedure, Interactive communication, and a cross reference to 
the various parameters for dynamic communication. 




Four Routines are used for basic program control. These control Routines 
are summarized on Table 3. 1.1.1. Routines START and PRINT are mandatory as 
they control Initialization, communication with boundary condition routines, 
timing, and output. Routines LOOP and CONVRG are optional depending on the 
configuration of the ECLSS model. 

Control and timing parameters are entered through interactive displays 
during active execution of START and PRINT. It may be desirable to dynamically 
communicate with the timing parameters during passive execution. This type of 
communication Is affected by entries In the user's MAIN. The following Control 
Routine descriptions Include Information on the Interactive displays as well 
as cross reference to the timing parameters. Examples of dynamic communication 
are Included In the sample model text. 


3. 1.1.1 
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Table 3. 1.1.1. Summary of Program Control Routines 
START Initialization control 

PRINT Boundary Condition Routine communication, timing 

control , and output control 

LOOP Timing control for first referenced ECLSS component 

In closed loop system 

CONVRG Convergence control In closed loop system 


3. 1.1.2 


3. 1.1.1 Routine START 


A flow diagram of the functions of Routine START Is given on Figure 
3.1. 1.1.1. 

A reference procedure Is 
CALL START 

for program execution Initialization. The call to Routine START Is the manda- 
tory first executable statement In the MAIN. 

Interactive communication with the Initialization parameters Is through 
console displays as shown on Figures 3. 1.1. 1.2 and 3. 1.1. 1.3. 

The cross reference to Initialization parameters Is shown on Table 3. 1.1. 1.1. 











RETURN 


Figure 3. 1.1. 1.1. Routine START (Concluded) 


3. 1.1. 1.3 
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3. 1.1. 1.4 
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5 INITIAL SYSTEM TEMP 621.100 DEG 





Figure 3.1. 1.1. 3. Typical Initialization Data Display 


3.1.1. 1.5 












3.1. 1.2 Routine PRINT 

A flow diagram of the functions of Routine PRINT Is given on Figure 
3.1. 1.2.1. 

The reference procedure Is 
CALL PRINT 

for program control. The call to PRINT Is the last Routine call In the timing 

loop. 

Basic Interactive communication Is through the console display shown on 

( 

Figure 3. 1.1. 2. 2. If the MAIN references component simulation routines which 
Imply consumables usage and/or orbital heating data Is required the display 
shown on Figure 3. 1.1. 2. 2 will be preceded by initialization and control data 
displays for these boundary condition functions. The Boundary Condition 
Routines are discussed In Section 3.2.1. Selection of the various output options 
will bring the various Output Routines and their associated control displays Into 
effect. The Output Routines are discussed In Section 3.2.2. 


3.1. 1-2.1 
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eRINT OPTION 


DISPLAY 

FLOG 

ITEM OPTION STATUS 

1 SELECT NODES O 

2 PLOTS 1 

3 SCHEMATIC 0 


************************************** 
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Figure 3. 1.1. 2. 2. Typical Print Option Display 
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| Q 3.1. 1.3 Routine LOOP 

Special timing control Is required for the nodal designation 1 associated 
with the first referenced ECLSS component In a closed loop system. Routine 
LOOP provides this control. 

The reference procedure 

CALL LOOP(I) 

which precedes the call to a component simulation routine with the nodal desig- 
nation I. 

There Is no further user communication with the functions or parameters of 
this Routine. The use of this routine Is mandatory for closed loop systems. 
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3. 1.1. 3.1 







3. 1.1. 4 Routine CONVRG 


Routines LOOP and CONVRG are used In conjunction for a closed loop system 
to control Iteration to a prescribed convergence tolerance. A typical applica- 
tion for Iteration about the node 1 follows: 


NN1 CONTINUE 
CALL LOOP(I) 

CALL PLATE(I.J) \ 

I 

( Typical component 
* ? simulation calls 

1 

CALL EVAP(M,I) / 

CALL C0NVRG(I,T0L,MAX,$NN2,$NN1) 

NN2 CONTINUE 


The call to CONVRG compares the Inlet fluid temperature at node I used In the 
first component call to the outlet fluid temperature at node I calculated by the 

last component call. If the two temperatures agree within the tolerance spec- 

• * * . m 

Ifled by the value of TOL, a return to statement number NN2 In the MAIN Is made. 
Whereas, If the agreement Is not In effect, the value of the temperature at 
node I for the first component call will be modified and the return will be to 
statement number NN1 In the MAIN. The parameter MAX Is a user selected maximum 
number of non-converging Iterations before CONVRG will terminate the execution. 

There Is no further user communication with the functions of Routine CONVRG. 
The use of Routine CONVRG Is optional. 
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3.1. 1-4.1 




3.1.2 
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3.1.2 Component Performance Routine 

This section provides Information on the simulation of various ECLSS 
components through application of the Library of Component Routines. This 
Library and a brief description of their function Is summarized on Table 
3.1. 2.1. 

Initial characteristics of the various components are entered or updated 
through an Interactive display during active execution. It Is often desirable 
to dynamically update or access various characteristics and variables during 
the passive execution. This latter updating Is accomplished by. entries In the 
user's MAIN which communicate with the subject routine variables. The following 
component routine descriptions Include Information on the Interactive dlsplpy 
as well as cross reference data for dynamic communication. Several simple 
dynamic updating examples are Included In the sample model. 



PLATE 

EVAP 

SPLIT 

MIX 

MOO 

RAO 

EXCH 

HEATER 

ATMO 

LIOH 

CONXG 


Table 3. 1.2.1. Summary of ECLSS Component 
Performance Routines 

Forced cooled Internal heat generating equipment 
Evaporator 

Branching of coolant flow 
Mixing of coolant legs Into a single junction 
Controlled proportioning of flow between two branches 
Radiator panel 

Counterflow and parallel flow heat exchanges 
Controlled fluid line heater 
Atmospheric compartment (Cabin) 

Carbon Dioxide removal with Expendable Lithium Hydroxide 
Condensing Heat Exchanger, condensing side 


CONXI 


Condensing Heat Exchanger, Interface (sink) side 


o 


3. 1.2.1 Routine PLATE 


The transient performance of forced cooled Internal heat generating equip 
ment Is simulated by Routine Plate. The routine Is an adaptive transfer 
function for the fluid flow segment from node 1 to J which processes the 
following equations. 

d T c 

( *>i • °i 1 • i e - < w) i AT i« • 

and 

<Vi 11 i * V • <*>i • 

where 




• V - < T c, 


L T 'i ’ Tj 


In 



The net loss of heat due to thermal coupling to m nodes defined by K(/) Is 
given by 


«c- £',<(/) < T ? ( 





where 


C 1K(/) * (ua) 1K(/) • 


and 


a ■ 1 


O 


for conduction or convective coupling, and. 


c « (/) ■ ctU mi) ’ 

and 

a - 4 

for radiation coupling of node 1 to K(£). 

The reference procedure Is 

CALL PLATE(I.J) 

to process the transfer function thiough node I to J. 

Interactive communication with the parameters and functions of the simulated 
component Is through a console display as Illustrated on Figure 3. 1.2, 1.1. An 
Integer value greater than zero assigned to Item 6, 7, and/or 8 will extend the 
communication upon exit from this display as follows: 

ITEM 6 Atmospheric Coolan t. The atmospheric parameters shown on Figure 

3. 1.2. 1.2 will be put through to node J. This feature Is used for 
continuity and mass conservation when PLATE Is used as a component 
simulation. 

ITEM 7 EPS Data Assignment . The Interactive display shown on Figure 

3. 1.2. 1.3 Is used to assign word numbers In the EPS data array to Q. . 

l 1 

(See Routine REPS.) 

ITEM 8 Thermal Coupling . The Interactive display shown on Figure 

3.1. 2. 1.4 Is used to thermally couple node I to m nodes defined by K (/). 

The 100 series nodes are used for convective coupling. The 200 series 
nodes are used for radiation coupling. To thermally couple node I to 
node 8 by conduction or convection assign 108 as the coupling code. 
Assignment of a 208 as the coupling will result In node I being coupled 


■, 


to node. 8 by radiation. Note the respective units of the coupling values 

V' . ' ’ “ "f* 

as displayed. 

- " -l 

The cross reference to Routine PLATE parameters Is shown on Table 3.1.2.111, 
Additional Information related to communication with the Thermal Coupling 
parameters Is glvep In the section on Routine COUPL. 


3.1. 2. 1.3 
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1 
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3 

4 
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DEG 
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ATMOSPHERIC COOLANT 

0 = NO 
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0 

INTEGER 

? 

EPS DATA ASSIGNMENT 
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0 
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8 
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0 

INTEGER 
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Figure 3. 1.2. 1.1. Typical Cold Plate Interactive Display 
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NOT EDITABLE 

PARTIAL PRESSURE OF WATER .154 PSIA 

PARTIAL PRESSURE OF NITROGEN 11.600 PSIA 

PARTIAL PRESSURE OF OXYGEN 3.100 PSIA 

PARTIAL PRESSURE OF CARBON .09? PSIA 

DIOXIDE 
ATMOSPHERIC 


PSIA 

PSIA 

PSIA 

PSIA 


PRESSURE 


1H.700 


H< >N sfc .\V \X' if: xV >*c XX >N :V \V \\* **if: X- * ^ xV xV ^ xV .{■ * \X >{< >}< ^ xV * ^ H< N v ^ 4* * 4* S< \V * {e X^ X^ * 


Figure 3. 1.2. 1.2. Typical Atmospheric Property Display 
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s$t * *• :* * jjr ;* *• ^ >;■: rf ■{< :fc -f ■*• *; # >N # : ';■: 


EPS DflTfi ASSIGNMENT 
FOR 

NODE NUMBER 6 

WORD 

ITEM NUMBER 

1 33 INTEGER 

2 3? INTEGER 

3 HD INTEGER 

* * * * * * * * * * * * * - * * ******** * * * * * 




«ijry 


Figure 3. 1.2. 1.3. Typical EPS Data Assignment Interactive Display 
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* ***y**«a^.. ..‘La-,. — — — .. * 7 ^ . . — 


c 


H/t sfc * % $ $ $ Jft tgt % «|e $ ifc ❖ * »{t * * *Jc sfc »> * H* ^ ate $$ ifc ifc ** ){« * >8< *8< $ ❖ a* »?< ♦ # 

THERMAL COUPLING DATA 
FOR NODE NO. 22 
COUPLED TO 1 NODES 
TYPE OF HEAT TRANSFER AND 
ITEM COUPLING NODE NUMBER CODE VALUE 
1 108 .800 
VALUES ARE: 

BTU/HR’FOR SERIES 100 COUPLING 
BTU/HR DEG**4 FOR SERIES 200 COUPLING 

Jfc ifrf ?ft >{< >fr: sfc >Jc >{; if; afc j{< >!« ^ >}• }$; ^ ^ ifc ^ ijt if s{: >$; ^ >}< *6t >f< sfc 



Figure 3. 1.2. 1.4. Typical Thermal Coupling Interactive Display 
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Table 3.1. 2.1.1. Dynamic Comnunl cation Cross Reference 
for Routine PLATE Parameters at Node I 


VARIABLE 

Component temperature 
Fluid Inlet temperature 
Thermal capacitance 
Heat transfer coefficient 
Coolant flow rate 
Internal heat generation 
Specific heat of gas mixture 
Partial pressure of HgO 

No 


Pressure of gas mixture 
Coupling values 1st ref. node 
2nd ref. node 



Tl 

(MC)i 

(UA)'j 

(WCp)i 


SYMBOL UNITS STORAGE I TYPE 


P( 

F< 

C( 

C( 

C( 

C( 

C( 

C( 

C( 

C( 

C( 

C( 
C(I,15 
C(I,16 


BTU/DEG 
BTU/HR DEG 
BTU/HR DEG 
BTU/HR 
BTU/LB DEG 
PSI 


;1K(1) 

: 1K(2) 


6th ref. node C^gj 
Number of coupled nodes m 

Coupling node number 1st ref. K(1) 

2nd ref. K(2) 


INTEGER 

INTEGER 

INTEGER 


C(I»20) 
IC 





6th ref. K( 6) INTEGER IC(I,10) R 



+ X « L present value 
X ■ K previous value 

*BTU/HR DEG for convection and conduction coupling (series 100) 
BTU/HR DEG 4 for radiation coupling (series 200). 
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3. 1.2.2 Routine EVAP 

The transient performance of an evaporator used as an ultimate heat sink 
Is simulated by Routine EVAP, The routine Is a transfer function for the fluid 
flow segment from node 1 to j which processes the following equations. 

(%)< < T j - V * <«>, 4T f. • 

where 



The media evaporated over the time span a to b Is calculated as 

^fg ^ ^p^l ^ 

and Is withdrawn from a source n as prescribed by the user. 

The simulation assumes the media to be at the saturation temperature T $ 
and does not account for sensible heat requirements to achieve this temperature. 
The reference procedure Is 

CALL EVAP(I,J) 

to process the transfer function through node I to J. 

Interactive communication with the parameters and functions of the simulated 
component Is through a console display as Illustrated on Figure 3. 1.2. 2.1. 
Reference to this Routine will automatically bring up the consumables source 
Initial loading display as Illustrated on Figure 3. 1.2. 2. 2. This display Is 
brought up Immediately prior to passive execution and Includes all referenced 
consumables sources as well as the source referenced through a particular call 
to EVAP. 


3.1. 2. 2.1 




c 

The cross reference to Routine EVAP parameters Is shown on Table 3.1. 2.2.1. 
Additional Information related to source assignment and Initial loading quantity 
for the evaporating media Is given In the section on Routine CONSUM. 
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cn-rcsi 


G 


NODE NUMBER 5 
EVAPORATOR 


TEM 


HEAT OF VAPORAZAT I ON 
OVERALL HEAT 
TRANSFER COEF. 
COOLANT FLOW RATE 
SATURATION TEMP 
INI TAL FLUID 
INLET TEMP 
CONSUMABLE 

1 = POT ABLE WATER 

2 = WATER 

3 - AMMONIA 

4 - OTHER 

TANK ASSIGNMENT 
INLET NODE NUMBER 
5 


10*0.000 BTU/LB 


3000.000 

1000.000 
495.000 


BTU/HR DEG 
BTU/HR DEG 
DEG 


.000 

1 


DEG 

INTEGER .LE. 


10 


6 INTEGER 
OUTLET NODE NUMBER 
& 


.LE. 20 




W 

>!' 


;^al MGe is 

■ ' w quality 


Figure 3. 1.2. 2.1. Typical Evaporator Interactive Display 
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ITEM 

SOURCE 

NO. 

CONSUMABLES 
SOURCES 
TYPE OF, 
CONSUMABLE 

UNIT 

INITIAL 

LOADING 

1 

1 

NITROGEN 

LBS 

.000 

2 

2 

OXYGEN 

LBS 

. 000 

3 

3 

LITHIUM HYDROXIDE 

LBS 

.000 

4 

4 

WATER 

LBS 

.000 

5 

5 

ELECTRIC FOWER 

WATT HRS 

.000 

6 

6 

POTABLE WATER 

LBS 

.000 


jJ( gj( ^ ^ tft }Q( ^ s{( )|( >|< )j/ }|< jj. ?}( ? ^t }f!){ i i )( n)t ?{ t){ , ^t '{'!f‘. )^( >f< ^ ^ ^ ^ ^ 




Figure 3. 1.2. 2. 2. Typical Consumables Source Display 
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Table 3. 1.2. 2.1. tynamlc Communication Cross Reference for 
tne Routine EVAP Parameters at Node I 


VARIABLE 

SYMBOL 

UNITS 

STORAGE 

TYPE 

COMMON 

Saturation temperature 

T , 

0E6 

p(i.* + > 

R 

MEAT 

Fluid Inlet temperature 

T i 

OEG 

F(I.X t ) 

R 

MEAT 

Heat of vaporatlon 

h fg 

BTU/LB 

C(I,1) 

R 

MEAT 

Heat transfer coefficient 

(«A) i 

BTU/HR DEG 

C(1 .2) 

R 

MEAT 

Coolant flow rate 

(WCp)i 

BTU/HR OEG 

C(1 .3) 

R 

MEAT 

Consumables type 

- 

INTEGER 

IC(I.3)* 

I 

MEAT 

Tank (source) number 

n 

INTEGER 

ICO .4) 

I 

MEAT 

Media consumed 

M 

LBS 

CONTOT(N) 

R 

CONS 


+ X * L present value 
X » K previous value 

*The consumables type Is used only for print out heading control and Is 
transferred to IAMCON(N) during Initialization. IAMC0N(N) Is stored In 
the common CONS. (See Routine CONSUM.) 
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3.1. 2.3 Routine SPLIT 


The splitting of flow at node 1 In n branches with a fixed flow proportion 
n J(K) at aach node 4(R) Is simulated by Routine SPLIT. The routine processes 
the equations 


" ®J(K) ***Vl 


T J(K) • T 1 


K ■ l,n 


The reference procedure Is 
CALL SPLIT(I t N) 


to split node I onto N branches. 

Interactive communication with the parameters and functions of the branching 
Is through a console display as shown on Figure 3.1. 2.3.1. An Integer value of 
unity assigned to Item 4 will bring up the Atmospheric Property display as shown 
on Figure 3. 1.2. 1.2. The atmospheric properties at node I will be put through 
to the branches j(K). 

The cross reference to Routine SPLIT parameters Is shown on Table 3. 1.2. 3.1. 



NODE NO. 12 
BRANCH 

SPLIT INTO 2 LEDS 
LED NODE FLOW 

ITEM NUMBER PROPORTION 

1 13 .500 

2 17 .500 

3 ATMOSPHERIC COOLANT 0 INTEGER 

0 = NO 

1 = YES 

H COOLANT FLOW RATE 769.090 BTU/HR DEG 

5 COOLANT INLET TEMP S20.317 DEG 


* * ♦ ♦ * * + + * # * M sf # + ■* * * * # # loft 


I 


f 

i 

f 



Figure 3. 1.2. 3.1. Typical Branching of Flow Interactive Display 
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Table 3. 1.2. 3.1. Dynamic Conmunl cation Cross Reference for 
the Routine SPLIT Parameters at Node I 


VARIABLE 

SYMBOL 

UNITS 

STORAGE 

Fluid Inlet temperature 

T i 

DEG 

F(I.X f ) 

Coolant flow rate 

<Vi 

BTU/HR DEG 

C(I. 3) 

1st branch node number 

JO) 

INTEGER 

ICO. 2) 

2nd branch node number 
• 

J(2) 

• 

INTEGER 

• 

ICO ,3) 
« 

• 

9th branch node number 

• 

J(9) 

• 

INTEGER 

• 

IC(IJO) 

1st branch flow proportion 

M J(1) 

FRACTION 

C0.12) 

2nd branch flow proportion 
• 

• 

H 4(2) 

• 

• 

FRACTION 

• 

• 

C(I, 13) 
• 

* 

9th branch flow proportion 

• 

M J(»> 

• 

FRACTION 

• 

CO, 20) 

Specific heat of gas mixture 


BTU/LB DEG 

C(I,5) 

Partial pressure of H 2 0 


PSI 

C(I,6) 

"2 


PS I 

CO. 7) 

°2 


PSI 

C0.8) 

CM 

8 

• 

PSI 

C( I ,9) 

Pressure of gas mixture 


PSI 

C(I,10) 


X ■ L present 
X » K previous 
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3. 1.2. 4 Routine MIX 


the mixing of flow from n legs defined by the nodes j(K) Into a single 
node 1 Is simulated by Routine MIX. The routine processes the equations 

“Vi • £ 

i£j“Vi(ic) T J(K) 

T 1 wp; — 

The reference procedure Is 
CALL MIX(I,N) 

to mix N legs Into the junction I . 

Interactive comnunl cation with the parameters and functions of the mixing 
Is through a console display as shown on Figure 3. 1.2. 4.1. An Integer value of 
unity assigned to Item 4 will result In processing of the atmospheric properties 
at the N nodes Into the node I. 

The cross reference to Routine MIX parameters Is shown on Table 3. 1.2. 4.1. 


C 


3. 1.2.4. 1 




c 


NODE NO. 5 
JUNCTION 
MIXING 3 NODES 
MIXED NODE 

ITEM NUMBER 

1 11 

2 16 

3 20 

H ATMOSPHERIC COOLANT 

0 = NO 

1 = YES 

* *: ifc ^ jfc sfc $ }{c i$c >fc 5$; :$< # ?M: # >}c >}c Jfc # Jfr 


UNIT 
INTEGER 
INTEGER 
INTEGER 
I INTEGER 


Figure 3. 1.2. 4.1. Typical Junction Interactive Display 
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Table 3. 1.2. 4.1. Dynamic Communication Cross Reference 
for Routine NIX at Node I 


VARIABLE 

SYMBOL 

UNITS 

STORAGE 

TYPE 

COMMON 

Fluid temperature at junction 

T i 

DEG 

F(I.X + ) 

R 

MEAT 

Coolant flow rate at junction 

( V.i 

BTU/HR DEG 

C(1 ,3) 

R 

MEAT 

1st mixing node number 

0(1) 

INTEGER 

MIXUM(H)* 

I 

MIXUP 

2nd mixing node number 

0(2) 

INTEGER 

MINUM 

(MO) 

I 

MIXUP 


nth mixing node number 

J(n) 

INTEGER 

MIXUM 

(M+N-l) 

I 

MIXUP 

Specific heat of gas mixture 


BTU/LB DEG 

CCI.5) 

R 

MEAT 

Partial pressure of H^O 


PSI 

C(I,6) 

R 

MEAT 

"2 


PS I 

C(I,7) 

R 

MEAT 

°2 


PSI 

C(1.8) 

R 

MEAT 

C0 2 


PSI 

Cd.9) 

R 

MEAT 

Pressure of gas mixture 


PSI 

0(1 .10) 

R 

MEAT 


r X ■ L present 
X « K previous 

*The value of M Is dynamically assigned during initialization and Is one 
greater than the total number of mixed nodes referenced prior to a 
particular call to MIX. M can be determined as the location of the value 
I In the IMIX array. That Is, IMIX(M) = I. IMIX Is In the common MIXUP. 
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3. 1.2. 5 Routine MOD 

The branching of flow at node 1 onto two legs defined by nodes j and n such 
that the flow proportioning f j and * n Is modulated to maintain a fixed tempera- 
ture T c at a control node m Is simulated by routine MOD. The routine processes 
equations similar to Routine SPLIT except that the flow proportioning Is 
dynamic* such that 

* f j * *< T » - v 

where aiDj Is a change In the flow proportioning factor from the previous 
calculation* g Is a proportional gains, and 



A maximum and minimum value for fj Is prescribed by the user. Note that a 
positive gains will favor node j if T n Is greater than T c (l.e., node j In this 
case Is assumed to be the cooling leg). 

The reference procedure Is 

CALL MOD(I) 

to modulate the branching of flow at node I. 

Interactive communication with the parameters and functions of the modulation 
Is through a console display as shown on Figure 3. 1.2. 5.1. An Integer value of 
unity assigned to Item 11 will bring up the atmospheric properties display for 
node I (see Figure 3. 1.2. 1.2) and put through the properties to the branches 
J and N. 

The cross reference to Routine MOD parameters Is shown on Table 3. 1.2. 5.1. 
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*****i|^****%********i|i)|^)|^i|c**iMl^il^*********i|c*********i|t*%**ifi 

NOOE NUMBER !0 
MODULATION VALUE 


O 


VALUE 

12 

11 

5 

505.000 

521.000 

.001 

l.i 


DESCRIPTION 
LEG 1 NODE NUMBER ' 

LEG 2 NODE NUMBER 
CONTROL NODE NUMBER 
CONTROL TEMP 

INITIAL TEMP AT CONTROL NODE 
PROPORTIONAL GAINS 
MAX HARD OVER 
MIN HARD OVER 
INITIAL TEMP AT MOD NODE 
COOLANT FLOW AT MOO NOOE 
ATMOSPHERIC COOLANT 

0 = NO, 

1 = YES 

* MUST BE DEFINED BEFORE YOU EXIT DISPLAY 

♦ * IQC * 9fc ♦ * * Ifc * # * 5fl« J#e 5je # afic $ ♦ j$t * ate * ♦ ♦ # aQc ifc # sfcHe ifc i|c # sjc * $ jfc % ?fc aQc aftc * * * * * jft afc * 


ITEM 

V* 

2 * 

3* 

4 

5 

6 

7 

8 
9 

10 

11 


518.891 


UNIT 

INTEGER 

INTEGER 

INTEGER 

DEG 

DEG 

FRACTION/DEG 
FRACTION 
FRACTION 
DEG 


1490. 48HBTU/HR. DEG. 
0 INTEGER 

\ 

\ 


i 


i 


i 



Figure 3. 1.2. 5.1. Typical Modulation Value Interactive Display 
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c 


Table 3.1. 2. 5.1. Dynamic Communication Cross Reference 
for Routine MOD at Node I 


VARIABLE 

SYMBOL 

UNITS 

STORAGE 

TYPE 

COMMON 

Leg 1 node number 

j 

INTEGER 

IC(I,4) 

I 

MEAT 

Leg 2 node number 

n 

INTEGER 

IC(I,5) 

I 

MEAT 

Control node 

m 

INTEGER 

IC(I,6) 

I 

MEAT 

Control temperature 

T c 

DEG 

C(I»15) 

R 

MEAT 

Temperature at control node 

T » 

DEG 

F(M,X + ) 

R 

MEAT 

Proportional gains 

9 

FRACTION 

C(I»16) 

R 

MEAT 

Maximum hard over 

MAX f j 

FRACTION 

C(I»17) 

R 

MEAT 

Minimum hard over 

MIN f. 

FRACTION 

C(I, 18) 

R 

MEAT 

Coolant temperature at MOD 
node 

T 1 

DEG 

F(I,X + ) 

R 

MEAT 

Coolant flow at MOD node 

<Vl 

BTU/HR DEG 

C(I,3) 

R 

MEAT 

Specific hear of gas mixture 

- 

BTU/LB DEG 

C(I,5) 

R 

MEAT 

Partial pressure of HgO 

- 

PSI 

C(I ,6) 

R 

MEAT 

n 2 

- 

PSI 

C(I »7) 

R 

MEAT 

°2 

- 

PSI 

C(I »8) 

R 

MEAT 

C 2 

- 

PSI 

C(I »9) 

R 

MEAT 

Pressure of gas mixture 

- 

PSI 

C(I,10) 

R 

MEAT 


+X * L present 
X * K previous 


C 

1 

* 
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3.1. 2.6 Routine RAD 

The transient performance of a radiator consisting of fluid flow paths 
attached to, or Integral with, a panel Is simulated by Routine RAD. Each fluid 
flow path Is paralled by a right and left section of the panel. The panel Is 
subject to Incident heat resulting from on-orbit operation. The Routine Is an 
adaptive transfer function for the fluid flow segment from node 1 to j which 
processes the following equations. 

d T c 

<*>1 ~3T * ' <«REJ * «c> ' < UA >1 4T f* • 
and 

<Vi (T j - v * (UA) 1 4T /« ’ 

where 



Qabs Is the absorbed heat, and Q c Is the net flow of heat resulting from 
thermal coupling. (The thermal coupling calculations are the same as those 
used for a cold plate. Section 3. 1.2.1.) 

The reference procedure Is 

CALL RAD(I,J) 

to process the transfer function through node I to J. Reference to this routine 


3.1. 2. 6.1 


will automatically bring up an Interactive display for definition of orbital 
Information Immediately prior to passive execution. (See Routine TRAJ.) 

Interactive communication with the parameters and functions Is through a 
console display as Illustrated on Figure 3. 1.2. 6.1. Positioning of the panel 
with respect to the spacecraft axes Is defined by the angle of Incident and 
dihedral angle as shown on Figure 3. 1.2. 6.2. The normal vector Is calculated 
Internally. An Integer value of unity assigned to Item 14 will bring up the 
Interactive node coupling display as for a cold plate. (See Section 3.1 .2.1 • 

Figure 3. 1.2. 1.4.) An Integer value greater than zero for Item 15 Indicates 
shadowing by the node number Indicated. The Interactive display for communication 
with the shadowing parameters shown on Figure 3. 1.2. 6. 3 Is brought up In this 
case. This display Is In communication with the shadowing routines. (See Routine 
SHAD.) 

The cross reference to Routine RAD parameters Is shown on Table 3. 1.2. 6.1. 
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NODE NUMBER 13 
RADIATOR PANEL 


ITEM 


VALUE 

UNIT 

• 

1 

THERMAL CAPACITANCE 

26.000 

BTU/DE6 

2 

OVERALL HEAT TRANSFER COEF. 

1500.000 

BTU/HR 0E6 

3 

COOLANT FLOW RATE 

384. S46 

BTU/HR DEG 

4 

SOLAR ABSORBI VITY 

. 100 

FRACTION 

6 

EMMISIVITY 

.900 

FRACTION 

6 

RIGHT FIN EFFECTIVENESS 

.800 

FRACTION 

7 

LEFT FIN EFFECTIVENESS 

.800 

FRACTION 

• 8 

RIGHT FIN AREA 

10.000 

SO FT 

9 

LEFT FIN AREA 

10. 000 

SO FT 

10 

ANGLE OF INC1DINCE 

.000 

RAD 

11 

DIHEDERIAL ANGLE 

.000 

RAD 

12 

INITIAL FIN TEMP 

621.000 

DEG 

13 

INITIAL COOLANT INLET TEMP. 

618.891 

DEG 

14 

NODE COUPLING 

0 

INTEGER 

IS 

NO = 0 
YES = 1 

SHADOW NODE NUMHER 

26 

INTEGER 


INLET NODE NUMBER OUTLET NODE NUMBER 


13 14 





i 

i 

1 

I 


L 


Figure 3. 1.2. 6.1. Typical Radiator Panel Interactive Display 
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*si<r»cn-ro4M»* 


SHADOWING DATA 
FOR NODE NO. 13 
SHADOWED BY 26 


ITEM 


SHADOW NODE AREA 


1.000 

$Q FT 

ANGLE OF INCIDENCE 


.000 

RAD 

DIHEDRAL ANGLE 


.000 

RAD 

STAND-OFF VECTOR DATA 


26 

INTEGER 

STANO-OFF DISTANCE 


2.000 

FT 

EQIIIV. STAND-OFF ANGLE OF 

INCIDENCE 

.000 

RAD 

EGUIV. STAND-OFF DIHEDRAL 

ANGLE 

.000 

RAO 


♦ * * * * * ** * ♦ ** * * * ** [* % % * ** ** if* * * * He * ** * * ** * ♦ * * * ♦ * * * * ♦ it* «Cc * * * 


Figure 3. 1.2. 6. 3. Typical Shadowing Parameter Interactive Display 


Table 3. 1.2. 6.1. Dynamic Communl cation Cross Reference 
for Routine RAD Parameters at Node I 


VARIABLE 

SYMBOL 

UNITS 

STORAGE 

TYPE 

COMMON 

Base temperature 

Tc 1 

DEG 

KM*) 

R 

MEAT 

Fluid Inlet temperature 

h 

DEG 

F(I.X) 

R 

MEAT 

Thermal capacitance 

(K), 

BTU/DEG 

C(I.l) 

R 

MEAT 

Heat transfer coefficient 

(UA) 1 

BTU/HR DEG 

C(1.2) 

R 

MEAT 

Coolant flow rate 

<Vi 

BTU/HR DEG 

C(I,3) 

R 

MEAT 

Incident heat 

C 

BTU/HR FT 2 

0(1.4) 

R 

MEAT 

Solar absorbability 

a 

FRACTION 

0(1,5) 

R 

MEAT 

Thermal emlsslvlty 

6 

FRACTION 

C(I .6) 

R 

MEAT 

Right fin effectlvlty 

r l 

FRACTION 

0(1.7) 

R 

MEAT 

Left fin effectlvlty 

r 2 

FRACTION 

C(I»8) 

R 

MEAT 

Right fin surface area 

A 1 

FT 2 

C(I»9) 

R 

MEAT 

Left fin surface area 

A 2 

FT 2 

C(I.IO) 

R 

MEAT 

Angle of Incidence 

a 

RAD 

C(I.ll) 

R 

MEAT 

Dihedral angle 

6 

RAD 

C(I»12) 

R 

MEAT 

Coupling value 1st ref node 

C 1K(1) 

* 

C(I»15) 

R 

MEAT 

2nd ref node 
« 

C 1K(2) 

• 

• 

* 

• 

C(I *16) 
• 

9 

R 

• 

MEAT 

• 

• 

6th ref node 

• 

C 1K(6) 

• ! 
* 

• 

C(I,20) 

• 

R 

• 

MEAT 

Number of coupled nodes 

M 

INTEGER 

IC(I,4) 

I 

MEAT 

Coupling node number 1st ref 

K(l) 

INTEGER 

IC( 1.5) 

I 

MEAT 

2nd ref 
• 

K(2) 

• 

INTEGER 

• 

IC(I *6) 
• 

I 

# 

MEAT 

• 

• 

6th ref 

* 

• 

K(6) 

• 

INTEGER 

• 

IC(I.IO) 

ft 

I 

• 

MEAT 

Shadowing node number 


INTEGER 

IC(I ,2) 

I 

MEAT 


+ X • L present value 
X ■ K previous value 


*BTU/HR DEGjfor convection and conduction coupling (series 100). 
BTU/HR DEG 4 for radiation coupling (series 200). 
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3. 1.2. 7 Routine EXCH 


( 

The exchange of heat In counter flow and parallel flow heat exchangers Is 
simulated by Routine EXCH. The routine Is a transfer function for the fluid 
flow segment from node 1 to J where the Interfacing fluid flow from m to n. 

The routine processes the equation for counter and parallel flow heat exchangers 
based of the methods of Reference 1. 

The reference procedure Is 

CALL EXCHO'J.M.N) 

to process the transfer function through node I to J based on the Interfacing 
condition at nodes M and N. A second call Is normally referenced In the Inter* 
facing coolant loop part of the model as 

( CALL EXCH(M,N,I ,J) 

which processes the transfer function through node M to N based on conditions at 
nodes I and J. The second call Is not mandatory. It Is used only as the model 
requires updating of the conditions at node N In the Interfacing loop. 

Interactive communication Is through a console display as Illustrated on 
Figure 3. 1.2. 7.1. Note that Information pertaining to the Interfacing side 
Is also Included. 

An Integer value of unity for Item 5 will cause the routine to put through 
the atmospheric properties from node I to J. This routine should not be used 
when the conditions at the Interface could result In condensation. (See Routines 
C0NX6 and CONXI . ) 

The cross reference to Routine EXCH parameters Is shown on Table 3. 1.2. 7.1. 

C 
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************.M*** + «**********.****«^*iM**** + 4>>****»**i|t*itii».*******it.******^ 

NODE NO. 21 
HEAT EXCHANGER 


ITEM 



VALUE 

UNIT 

• 


*+* CALLING 

side *** 



• 

1 

HEAT TRANSFER COEF. 


300. 000 

BTU/HR DEG 

• 

2 

COOLANT FLOW RATE 


100. 000 

BTU/HR DEG 

• 

3 

FLUID INLET TEMPERATURE 


52 1.000 

DEG 

• 

4 

TYPE 


0 

INTEGER 

• 

5 

COUNTERFLOW s 0 
PARALLEL FLOW = 1 
ATMOSPHERIC COOLANT 


0 

INTEGER 

• 


NO r 0 
YES r l 

*<♦ INTERFACE 

SIDE *** 



• 

S 

HEAT TRANSFER COEF. 


3000. 000 

BTU/HR DEG 

• 

7 

COOLANT FLOW PATE 


14$»0. 484 

BTU/HR DEG 

• 

e 

FLUID INLET TEMPERATURE 


518. 751 

DEG 


4 + <' 


:**>:>* 


* * * * * * * * * *< \*< * 4 * f * K K *i* 4 if' A »! * A- A > ;! < * * 


A- -I A * *fr A' A A A< 1 K- *: 4 »{ *« ♦ »? A .J rf. 


* ; A 


r 1gur« 3. 1.2. 7.1. Typical Counter or Parallel Flow 

Heat Exchanger Interactive Display 
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Table 3. 1.2. 7.1. Dynamic Communication Cross Reference for Routine 
EXCH Parameters at Node I Interfaced to Node M 


VARIABLE 

SYMBOL 

UNITS 

STORAGE 

TYPE 

COMMON 

Calling Side 






Fluid Inlet temperature 

T i 

nrp 

.. U 

F(I»X + ) 

R 

MEAT 

Heat transfer coefficient 

(UA) i * 

BTU/HR DEG 

C(I,2) 

R 

MEAT 

Coolant flow rate 

»Vi 

BTU/HR DEG 

C( I ,3) 

R 

MEAT 

Specific heat of mixture 

- 

BTU/LB DEG 

C(I»5) 

R 

MEAT 

Partial pressure of H 2 0 

- 

PSI 

C(I »6) 

R 

MEAT 

n 2 

- 

PSI 

C( I ,7) 

R 

MEAT 

°2 

- 

PSI 

C(I ,8) 

R 

MEAT 

C0 2 

- 

PSI 

C( I ,9 ) 

R 

MEAT 

Pressure of mixture 

- 

PSI 

C(I,10) 

R 

MEAT 

Interface Side 






Fluid inlet temperature 

T m 

DEG 

F(M,X + ) 

R 

MEAT 

Heat transfer ceofficient 

<«>■* 

BTU/HR DEG 

C(M,2) 

R 

MEAT 

Coolant flow rate 


BTU/HR DEG 

C(M,3) 

R 

MEAT 






+ X = L present value 
X = K previous value 

h * ® 

The Routine processes an overall heat transfer coefficient 

< UA >0 ■ 1 7 - r- 

TWi TuAlT 
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3. 1.2. 8 Routine HEATER 



IT 


The performance of an In-line fluid heater, whose power Is Qp^ switched 
In response to the fluid temperature at a control node n and a control tempera- 
ture of T c . The routine Is an adaptive transfer function for the fluid flow 

segment 1 to j and processes equations as for Routine PLATE except that Q. , Is 

*1 

either Q pwR or zero depending on the on/off configuration, respectively. 

The on/ off configuration is control to a dead-band aT about the control 

temperature T r . The heater Is switched on when T n Is less than T - aT and does 
c n c 

not go off until T n Is greater than T c + aT. The energy In watt hours consumed 
over the time span a to b Is calculated as 

r b 

E * 3.4130 / Q. dt. 

-'a 

and Is drawn from a source lx as prescribed by the user. 

The reference procedure Is 

CALL HEATER(I.J) 

to process the transfer function through node I to J. 

Interactive communication with the parameters and functions of the In-line 
heater Is through a console display as Illustrated on Floure 3. 1.2.8. 1. rnteoer 
values of unity for Items 10 and/or 11 will put through the atmospheric prop- 
erties and/or set up for thermal coupling data entry, respectively. (See 
Routine PLATE). Reference to this Routine will automatically bring up the 
consumables source Initial loading display as Illustrated on Figure 3. 1.2. 2. 2. 
This display Is brought up Immediately prior to passive execution and Includes 
all referenced consumables sources as well as the source referenced through a 
particular call to HEATER. 


3.1.2.P.1 


G 


The cross reference to Routine HEATER parameters Is shown on Table 3.I.2.8.I. 
Additional Information related to source assignment for the energy Is given In 
the section on Routine CONSUM. 




f 
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*******************************************+*****+********************* 

NODE NUMBER H 


a> 


HEATER 



ITEM 


VALUE 

UNIT 


1 

THERMAL CAPACITY 

200. 000 

BTU/DE6 


2 

OVEALL HEAT TRANSFER COEF. 

1500.000 

6TU/HR DEG 


3 

COOLANT FLON RATE 

SOI. 000 

BTU/HR DEG 


H 

INITIAL COMPONENT TEMPERATURE 

52 1.000 

OEG. 


5 

INITIAL INLET TEMPERATURE 

500.830 

DEG. 


6 

HEATER POMER 

450.000 

BTU/HR DEG 


? 

CONTROL NODE NUMBER 

2 

INTEGER 


8 

CONTROL TEMPERATURE 

532. OOC 

DEG. 


9 

DEAD BAND 

2.000 

DEG. 


10 

INITIAL TEMP AT CONTROL NODE 

510. 185 

DEG 


11 

ATMOSPHERIC COOLANT 

1 

INTEGER 



0 = NO 

1 = YES 




12 

THERMAL COUPLING 

0 

INTEGER 



0 = NO 

1 = YES 




13 

POWER SOURCE 

5 

INTEGER 



INLET NODE NO. 

OUTLET NODE 

NO. 


4 1 

***.; !j H :) ■e'4-***************«:******************>M'***«:**>Si*i»i************ 


Figure 3. 1.2. 8.1. Typical In-line Heater Interactive Display 
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Table 3. 1.2. 8.1. Dynamic Communl cation Cross Reference for 
Routine HEATER Parameters at Node I 


VARIABLE 


Component temperature 
Fluid Inlet temperature 
Thermal Capacitance 
Heat transfer coefficient 
Coolant flow rate 
Power Applied 
Control node number 
Control temperature 
Deadband 

Temperature at control node 
Heater power (ON) 

Energy source 
Energy 

Specific heat of mixture 
Partial pressure of HgO 


SYMBOL UNITS STORAGE TYPE COMMON 


Pressure of mixture 
Coupling value 1st ref. node 
2nd ref. node 


'1 

(MC) i 

(UA) i 

(WCp)i 


DEG 

DEG 

BTU/DEG 
BTU/HR DEG 
BTU/HR DEG 
BTU/HR 
INTEGER 
DEG 
DEG 
DEG 

BTU/HR 
INTEGER 
WATT HRS 
BTU/LB DEG 
PSI 


;ik(i) 

'1K(2) 


6th ref. node * 

Number of coupled nodes M INTEGER 

Coupling node number 1st ref. K(1) INTEGER 

2nd ref. K(2) INTEGER 


6th ref. 


INTEGER 



C(I»4) 
IC(I,2) 
C(I,11) 
C(I,13) 
F(N,X + ) 
C(I»12) 

IC(I ,3) 
CONTOT(IX) R 
C(I ,5) 

C( 1 ,6) 

C(I »7) 

C(I,8) 

C(I,9) 

C( I ,10) 
C(I,15) 
C(I»16) 


C(I .20) 
IC(I. 
IC(I. 

icci.i 



1 » L present value 
X'K previous value 

*BTU/HR DEG.for convection and conduction coupling 
BTU/HR DEG 4 for radiation coupling (series 200). 

3.1. 2. 8.4 


IC(I, 


(series 100). 




















3. 1.2. 9 Routine ATMO 


Processing of a gas stream through atmospheric compartment subject to heat 

addition, moisture generation, carbon dioxide veneration, oxypen consumption, and 

external leakage is simulated by Routine ATMO. The heat addition results from 

equipment heat generation, Q T , and the sensible metabolic heat of up to six 

l 1 

crew members. Crew sensible heat, water addition, carbon dioxide addition, and 
oxygen consumption are proportional to the user specified metabolic load of the 
occupants. Oxygen and Nitrogen external leakage make up Is supplied from 
sources n 1 and n ? specified by the user. The routine Is an adaptive transfer 
function for the gas flow segment from the inlet 1 to the exit j and performs 
a mass balance on the atmospheric constituents. The atmospheric conditions at 
node j represent the compartment conditions. 

The reference procedure is 

CALL ATMO (1,0) 


to process the transfer function from node I to J. The first component simulation 
call In an atmospheric loop should reference this routine. 

Interactive communication with the parameters and functions of the atmospher> 
1c compartment simulation is through console displays as shown on Figures 
3. 1.2. 9.1 and 3. 1.2. 9. 2. Reference to this routine will automatically bring up 
the consumables source initial loading display as Illustrated on Figure 3. 1.2. 2. 2. 
This display Is brought up Immediately prior to passive execution and Includes 
all referenced consumables sources as well as the source referenced through a 
particular call to ATMO. 

The cross reference to Routine ATMO parameters Is shown on Table 3. 1.2.9. 1 
Additional information related to source assignment for the oxygen and nitrogen 
is given in the section on Routine CONSUM. 
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NODE 1 

ATMOSPHERIC COMPARTMENT 


ITEM 


VALUE 

UNIT 


l 

COMPARTMENT VOLUME 

1000.000 

CUBIC FT 


2 

LEAKAGE RATE 

2.000 

LB/HR 


3 

COOLANT FLOW RATE 

SOI. 000 

BTU/HR DEG 


4 

HEAT LOAD 

2001.000 

BTU/HR 


S 

SPECIFIC HEAT OF GAS 

.210 

6TU/LB DEG 


6 

PARTIAL PRESSURE OF WATER 

.130 

PSI 


7 

PARTIAL PRESSURE OF NITROGEN 

11. *00 

PSI 


8 

PARTIAL PRESSURE OF OXYGEN 

3. 100 

PSI 


9 

PARTIAL PRESSURE OF CARBON DIOXIDE 

.093 

PSI 


10 

TOTAL PRESSURE 

14.700 

PSI 


11 

NITROGEN TANK 

1 

« 

INTEGER 


12 

OXYGEN TANK 

2 

INTEGER 


13 

INLET GAS TEMPERATURE 

SOS . 100 

DEG 



Figure 3. 1.2. 9.1. Typical Atmospheric Compartment Interactive Display 
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#****#5fc*#:*>M<**:****>M:**>M >i ** * .^ ********:**** >M * t* ******* >t< *#*****$*****$* *; 

CREW MEMBER 
METABOLIC RATES 


NODE NO. 1 


ITEM 

METABOLIC 

UNIT 


RATE 


1 

. 000 

BTU/HR 

O 

4 

. 000 

BTU/HR 

3 

.000 

BTU/HR 

4 

BOO. 000 

BTU/HR 

5 

. 000 

BTU/HR 

6 

. 000 

BTU/HR 

* *{« >;< * & ft >$< >{c *c ft # #: afc ft. >* >{ 

$: K :*•#**# * *************:** * :>}:*:**;**#* ^ ft >v $ ft * jjt * : * $ }: sfrr 


Figure 3. 1.2. 9. 2. Typical Crew Metabolic Data Interactive Display 
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Table 3, 1,2,9. 1 Dynamic Communication Cross Reference 

Parameters for ATMO Parameters at Note I 


VARIABLE 

SYMBOL 

UNITS 

STORAGE 

TYPE 

COMMON 

Gas Inlet temperature 

T i 

DEG 

Ft I. It*) 

R 

MEAT 

Compartment volume 

v i 

FT 3 

e(i.i) 

R 

MEAT 

Leakage rate 

M e1 

LB/HR 

C(I,2) 

R 

MEAT 

Gas flow rate 

<*p>1 

LB/HR DEG 

C(I,3) 

R 

MEAT 

Equipment heat load 


BTU/HR 

C(I,4) 

R 

MEAT 

Specific heat of gas mixture 

- 

BTU/LB DEG 

C(I,5) 

R 

MEAT 

Partial pressure of HgO 

- 

PSI 

C(1 ,6) 

R 

MEAT 

"2 

- 

PSI 

C(I,7) 

R 

MEAT 

°2 

- 

PSI 

C(I ,8) 

R 

MEAT 

co 2 

- 

PSI 

C(I,9) 

R 

MEAT 

Pressure of gas mixture 

- 

PSI 

C(,10) 

R 

MEAT 

Metabolic rate 1st member 

- 

BTU/HR 

C(I,15) 

R 

MEAT 

2nd member 
• 

• 

BTU/HR 

• 

C (I ,16) 
• 

R 

• 

MEAT 

• 

• 

6th member 

• 

• 

BTU/HR 

• 

C(1.20) 

• 

R 

• 

MEAT 

Source of N 2 

"l 

INTEGER 

IC(1 .5) 

I 

MEAT 

Same for 0 2 

n 2 

INTEGER 

ICO, 6) 

I 

MEAT 

N 2 consumed 

^N2 

LBS 

CONTOT(Nl) R 

MEAT 

0 2 consumed 

M 02 

LBS 

C0NT0T(N2) R 

MEAT 


+ X « L present value 
X ■ K previous value 




t 


i 



3.1.2.10 Routine LI0H 

The removal of Carbon Dioxide from a gas stream by a Lithium Hydroxide 
canister Is simulated by Routine LIOH. The routine Is an adaptive transfer 
function for the gas stream segment from node 1 to j which processes the follow- 
ing equations. 

p co 2 j 8 P C0 2 1 

where 4 >, an efficiency factor Is given by 


1 - e 


-a 


1 - e’ a + e 


.ae-a 


here 

K M, p, 

a * "Tf » 

W 1 

where K Is an empirically determined reaction rate (Reference 1) given by 

K • 1100. - 700. (1 - e 10 * 6 ) 0. < e < .8 

K * 400. 6 > .8 


and 



c co 2 W 1 


K, 


1 


dt 


where C CQ Is the Carbon Dioxide concentration. Heat and moisture Is added to 
the stream by the reaction. 

The amount of Lithium Hydroxide consumed Is equal to the number of canister 
changes time the mass of the canisters. A canister Is automatically changed 
when 





j 


b' 


i 
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P C0 2 j > (P C0 2 jWx- 

The Lithium Hydroxide Is withdrawn from a source n specified by the user. The 
reference procedure Is 

CALL LIOH(I.J) 

to process the transfer function through node I to 0. 

Interactive communication with the parameters and functions of the canister 
Is through a console display as shown on Figure 3.1.2.10.1 Reference to this 
routine will automatically bring up the consumables source Initial loading 
display as Illustrated on Figure 3. 1.2. 2. 2. This display Is brought up Inmedl- 
ately prior to passive execution and Includes all referenced consumables sources 
as well as the source referenced through a particular call to LIOH. The routine 
automatically processes the remaining atomospherlc properties. 

The cross reference to Routine LIOH parameters Is shown on Table 3.1.2.10.1. 
Additional Information related to source assignment for the Lithium Hydroxide 
Is given In the section on Routine CONSUM. 


C 
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NODE NO 


2 


* 


LITHIUM HYDROXIDE 
CANISTER 


ITEM 


VALUE 

UNITS 

1 

CANISTER MASS 

1.000 

LBS 

2 

CAS FLON RATE 

SOI. 000 

BTU/HR DEG 

3 

CANISTER PRESSURE CHANGE 

. ISO 

PSI 

H 

SPECIFIC HEAT OF GAS 

.210 

BTU/LB DEG 

S 

PARTIAL PRESSURE OF MATER 

. 130 

PSI 

6 

PARTIAL PRESSURE OF NITROGEN 

11.600 

PSI 

7 

OXYGEN 

3. 100 

PSI 

8 

CARBON 

. 093 

PSI 


DIOXIDE 



9 

TOTAL PRESSURE 

1H. 700 

PSI 

' 10 

INITIAL ABSORBED QUANTITY 

.000 

FRACTION 

*11 

CANISTER SOURCE 

3 

INTEGER 

' 12 

INLET GAS TEMPERATURE 

S10. 16S 

DEG 




Figure 3.1.2.10.1. Typical Lithium Hydroxide 

Canister Interactive Display 
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Table 3.1.2.10.1. Dynamic Cross Reference to 

Routine LI0H Parameters at Node I 


VARIABLE 

SYMBOL 

UNITS 

STORAGE 

TYPE 

COMMON 

Gas Inlet temperature 

T i 

OEG 

FO.X*) 

R 

MEAT 

Canister mass 

"i 

DEG 

Cll.11 

R 

MEAT 

Gas flow rate 

<v« 

BTU/HR DEG 

C(I ,3) 

R 

MEAT 

Canister change pressure (P { 

CoJmAX 

PSI 

CC.11) 

R 

MEAT 

Heat added 

4 

BTU/HR 

CO. 4) 

R 

MEAT 

Specific heat of mixture 

- 

BTU/LB DEG 

C(I »5) 

R 

MEAT 

Partial pressure of H 2 0 

- 

PSI 

C( I ,6) 

R 

MEAT 

n 2 

- 

PSI 

C(I »7) 

R 

MEAT 

°2 

m 

PSI 

C(I,8) 

R 

MEAT 

co 2 

p co ? 

PSI 

C(I#9) 

R 

MEAT 

Pressure of gas mixture 


PSI 

C(I,10) 

R 

MEAT 

Lithium Hydroxide source 

n 

INTEGER 

IC(I ,5) 

I 

MEAT 

Consumed Lit Hum Hydroxide 


LBS 

CONTOT(N) 

R 

CONS 


+ X ■ L present value 
X ■ K previous value 


3.1.2.19.A 





3.1.2.11 Routine C0NX6 


The performance of a counterflow condensing heat exchanger Is simulated by 
Routines CONXG and CONXI. Routine CONXG, discussed In this section, processes 
the condensing side of the heat exchanger. Routine CONXI, discussed In the 
next section, processes the Interface (or sink) side of the heat exchanger. 

The routine discussed In this section Is a transfer function for the gas stream 
segment 1 to j Interfacing with a counterflow segment m to n. Processing 1$ as 
follows and Is Illustrated on Figure 3.1.2.11.1. 

The gas and water vapor mixture enter at a temperature t^ and partial 

pressure P H The mixture experiences only sensible cooling until the dew 

point Is reached. The remaining portion of the heat exchanger dehumlfles and 

cools the mixture to tj and P H qj at the exit. The sensible cooling portion of 

the heat exchanger Is referred to as the "dry" section with a heat transfer 

coefficient of UA n If the entire heat exchanger were dry. The dehumldlfylng 
U 1 

section Is referred to as the "wet" section with a heat transfer coefficient 
of UAy If the entire heat exchanger were wet. Internal calculation proportion 
the dry and wet sections and the applicable value of the heat transfer coeffi- 
cients. The condensation Is stored In a source (tank) n^ as prescribed by the 
user. 

The reference procedure Is 
CALL CONXG ( I ,J M,N) 


to process the condensing side from node I through node J with nodes M and N as 
the respective Inlet and outlet of the Interfacing sink fluid. A call to CONXI 
Is mandatory In the Interfacing cov^ant loop part of the model. 

Interactive commuiil cation with the parameters and functions of the condensing 
side of the heat exchanger Is through a console display as shown on Figure 


3.1.2.11.1 


3.1.2.11.2. Reference to this routine will automatically bring up the consumables 
source Initial loading display as Illustrated on Figure 3. 1.2.2. 2. This display 
Is brought up Immediately prior to passive execution and Includes all reference 
consumables sources as well as the source referenced through a particular call 
to CONXG. The routine automatically processes the remaining atmospheric 
properties. 

The cross reference to Routine CONXG parameters Is shown on Table 3.1.2.11.1. 
Additional Information related to storage assignment for the condensation Is given 
In the section on Routine CONSUM. 
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wmmmmmm 






NODE NO. 3 

CONDENSING HEAT EXCHANGER 
ATMOSPHERIC SIDE 


‘ITEM 

1 

•> 

3 

4 



VALUE 

*** CALLING SIDE *** 


CONDENSING HEAT TRANSFER COEF. 
DRY HEAT TRANSFER COEF. 

COOLANT FLOW RATE 
FLUID INLET TEMP 
CONDENSATE TANK NO. 


2000. 000 
1000.000 
SOI . 000 
S 1 2. 278 
4 


*** INTERFACE SIDE *** 

COOLANT FLOW PATE 1490.484 

FLUID INLET TEMP 498.000 

CALLING SIDE NODES IN 3 OUT 4 

INTERFACE SIDE NODES IN 7 OUT 8 


UNIT 


BTU/HR DEG 
BTU/HR DEG 
BTU/HR DEG 
DEG 

INTEGER 


BTU/HR DEG 
DEG 


ft********* * ❖ * * * * * ❖ * * * * * * * >N * * X * * * * * * * * * * XV ****** * * K * * * xV * * ***** * * * * * * l 


■> K 


Figure 3.1.2.11.2. Typical Condensing Heat Exchanger 

(Condensing Side) Interactive Display 
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Table 3.1.2.11.1. Dynamic Cross Reference to Routine 

CONXG at Node i Interfaced to Node M 


VARIABLE 

SYMBOL 

UNITS 

STORAGE 

TYPE 

COMMON 

Gas Inlet temperature 

T i 

DEG 

F(I,X + ) 

R 

MEAT 

Heat transfer coefficient, 
wet 

u \ 

BTU/HR DEG 

C(I.l) 

R 

MEAT 

Heat transfer coefficient, 
dry 


BTU/HR DEG 

C( I ,2) 

R 

MEAT 

Gas flow rate 

<Vi 

BTU/HR DEG 

C( I ,3) 

R 

MEAT 

Heat rejected 

- 

BTU/HR 

C(H ,4) 

R 

MEAT 

Specific heat of mixture 

- 

BRU/LB DEG 

C( I ,5) 

R 

MEAT 

Partial pressure of H 2 0 

P H J> 01 

PSI 

C( I ,6) 

R 

MEAT 

N 2 

c 

PS I 

C(I ,7) 

R 

MEAT 

°2 

- 

PSI 

C(I ,8) 

R 

MEAT 

co 2 

- 

PSI 

C( 1 ,9) 

R 

MEAT 

Pressure of mixture 

- 

PSI 

C(I.IO) 

R 

MEAT 

Condensate tank number 

n l 

INTEGER 

IC(I ,3) 

I 

MEAT 

Condensate stored 

- 

LBS 

CONTOT(NI) R 

CONS 


X « L present value 
X * K previous value 


3.1.2.11.5 
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3.1.2.12 Routine CONXI 

Routine CONXI Is & mandatory companion routine to C0NX6 of the preceding 
section for simulation of a counterflow condensing heat exchanger. The 
routine discussed in this section Is a transfer function for the Interfacing 
side of a condensing heat exchanger fluid segment m to n with counterflow on 
the condensing segment from 1 to J. The routine processes the following equation 

^n " \i + ^n/^ WC p^m 

where Q m has been assigned by the processing of Routine CONXG. 

Interactive communication with the parameters and functions of the Inter- 
face side of a condensing heat exchanger Is through a console display as shown 
on Figure 3.1.2.12.1. Entry of an Integer value of unity for item 3 will 
put through the atmospheric properties of node m to node n. 

The cross reference to the parameters for interfacing side of a condensing 
heat exchanger are shown on Table 3.1.2.12.1. 
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0 


ftftftft ft ftcftftftfttftft ftfttft ftft ft'fttftftftftftftftft ftft ft fttft ftcft ftftftftfttftft ftftftftftft ftftftftftftftftftftftftftftftftftftftftftftfttft 

NODE MO. 7 

CONDENSING HEAT EXCHANGER 
INTERFACE SIDE 


ITEM 


1 

2 


«■** CALLING SIDE *** 


COOLANT FLOW RATE 
FLUID INLET TEMP. 
ATMOSPERIC COOLANT 
NO = 0 
YES = 1 


VALUE 


l»r«0.48H 
4'«S. 715 
0 


UNIT 


BTU^HR DEG 
DEG 

INTEGER 


*** ATMOSPHERIC SIDE **+ 

4 COOLANT FLOW PATE Mil. OOO 

5 FLUID INLET TEMP • • 1 ;* .276= 

CALLING SIDE NODES IN 7 ntJT 8 

ATMOSPHERIC SIDE NODES IN 3 "UT 4 

ft ft ft ft ft ft: ft ft ft ft: ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft: ft ft ft: ft ft ft ft ft: ft ft :«* ft ft ft ft :}: ft ft ft ft ft: ft ft ft ft ft ft: ft ft: ft: ft ftftt ft ft 


BTU/HR DEG 
DEG 


j 


i 


8 


V 


o 

Figure 3.1.2.12.1. Typical Condensing Heat Exchanger 

(Interface Side) Interactive Display 
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Table 3.1.2.12.1. Dynamic Cross Reference to 

Routine CONXI at Node M 


VARIABLE 

SYMBOL 

UNITS 

STORAGE 

TYPE 

COMMON 

Fluid Inlet temperature 

T m 

m 

DEG 

FC.X*) 

R 

MEAT 

Coolant flow rate 

<v« 

BTU/HR DEG 

C(M,3) 

R 

MEAT 

Heat absorbed 


BTU/HR 

C(M,4)* 

R 

MEAT 

Specific heat of mixture 

- 

BTU/LB DEG 

C(M,5) 

R 

MEAT 

Partial pressure of H 2 0 

- 

PSI 

C(M,6) 

R 

MEAT 

N 2 

- 

PSI 

C(M,7) 

R 

MEAT 

°2 

- 

PSI 

C(M,8) 

R 

MEAT 

co 2 

• 

PSI 

C(M,9) 

R 

MEAT 

Pressure of mixture 

- 

PSI 

C(M,10) 

R 

MEAT 


+ X - L present value 
X * K previous value 

*Va1ue assigned by companion call to CONXG. 




Q 3.1.3 Input Utility Routine 

Two routines are available for tabular data look-up and Interpolation. The 
routines TABLE and STEP as summarized on Table 3. 1.3.1 differ only In the manner 
of Interpolation. 

The Input Utility Routine descriptions which follow Include reference 
procedure and Interactive communication In the table data set up. 


c 


( 


Table 3. 1.3.1. Summary of Input Utility Routines I 

‘I 

) 

TABLE Linear Interpolation of dependent variable 

In tabular data. 

STEP Interpolates tabular data as step function. 


u ■ 

P - 

t: 

P 

S- 


o 


3.1. 3.2 
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3.1. 3.1. Routine TABLE 


Linear Interpolation of Input tabular data Is performed by Routine TABLE. 
The reference procedure Is 

CALL TABLE (NTAB, XXI ,YY1 ) 

where NTAB Is the table number. The Routine returns the linear Interpolated 
value of YY1 for the prescribed Independent variable XXI from the tabular arrays 
XX and YY. 

Interactive communication Is through a console display as shown on Figure 
3.1.3. 1.1. 


( 



3.1. 3.1.1 
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************************** 
TABLE 2 
XX YY 

1 .000 200.000 

2 2.S00 400. 000 

*###*##*####****#*###*###* 


Figure 3. 1.3. 1.1. Typical Interactive Tabular Data Display 


3.1. 3. 1.2 


34*3.2 Routine STEP 






Interpolation of Input step function tabular data Is performed by Routine 

STEP. 

The reference procedure Is 
CALL STEP(NTAB,XX1 ,YY1 ) . 

The functions and Interactive communl cation with the routine Is the same as for 
Routine TABLE except for the method of Interpolation. 


i 

i 
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i 
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3.2 UNREFERENCED ROUTINES 



( 



This section presents a description of the functions and parameters of those 
library routines which are automatically brought Into execution as a result of 
use of and/or specific option selection with respect to the referenced routines 
of Section 3.1. Only those library routines with which the user has Interactive 
communication and/or may desire dynamic communication are Included. These types 
of routines Include Boundary Condition Routines and Output Routines. 


3 . 2.1 






3.2.1 Boundary Condition Routines 

Five routines are directly associated with development of Boundary 
Conditions. The function of these routines are summarized on Table 3.2.1. 1. 


5 >’ 


I 

l 




Table 3. 2. 1.1. Summary of Boundary Condition Control Routines 

' y 

COUPL Provides for Initial cross coupling of thermally 

connected nodes. 

REPS Provides data array for electrical power assignment. 

CONSUM Provides for Integration of expended or generated 
media for various source assignments. 

TRAJ Performs Incident (orbital) heating calculations. 

SHAD Assembles absorbed (orbital) heat fcr a panel. 



3.2. 1.2 
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I 

3.2*1. 1 Routine COUPL 

All thermo 1 coupling date Is performed Internal to the argumented component 
performance routines* Initial coupling data Is loaded through the Interactive 
displays associated with these routines. Routine COUPL performs the cross 
coupling during this Initial loading. If the user specifies node I Is coupled 
to node J, Routine COUPL automatically couples node J to I with the same coupling 
value. This function Is performed only during Initial loading. Accordingly, 
any dynamic communication with coupling values and/or type of coupling requires 
that the user search out the location and change the coupling data for both 
node I and J. 




3.2.1. 1.1 


3.2. 1.2 Routine REPS 

* 

Electric power date assignment Is controlled through the Initial Interactive 
communication with the argument component performance routine reference. Routine 

WPS provides the Interface of the data assignment to the electrical power data 

\ 

to be read from tape or assigned through an Interface to another program. The 
communication Is trhough the dimensioned variables NEPS (100) and DEPS(IOO) In 
the COMMON/EPS/. The resident value P 1 of electrical power (BTU/HR) In DEPS(I) 

Is assigned to node J by setting 
NEPS(I) ■ J. 

The argumented component performance routine at node J searches the NEPS array 
for the Integer J and sums all corresponding values of OEPS into C(J,4). Dynamic 
reassignment Is accompllsjed by control of the NEPS array. 

The current version of Routine REPS assigns specific Shuttle Prfciter 
electrical power heating value to the NEPS array words shown on the menu Illus- 
trated on Figure 3.2. 1.2.1. The routine reads a dictionary of active components 
from Unit 10 and the power timeline from Unit 11. 


3.2. 1.2.1 
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MPL3 

> 

IF YOU WONT TO DISPLAY MENU ENTER 1, OTHERWISE BLANK 

>1 

STS 

HEAT LOAD MENU 


WORD NO. 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

PAUSE FOR 

> 


DESCRIPTION 

AVIONICS AIR BAY 1 

AVIONICS AIR BAY 2 

AVIONICS AIR BAY 3 

AVIONICS FAN 1A/18 BAY =1 

AVIONICS FAN 2A/2B BAY =2 

AVIONICS FAN 3A/3B BOY =3 

CABIN AIR COOLED 
CABIN FAN 

COLDPLATE FREON DFI MID-BDY CONTAINER =1 
COLDF'LATE FREON DFI MID-BDY CONTAINER =2 
COLDPLATE FREON DFI MID-BDY CONTAINER x3 
TBD 

COLDPLATE WATER DFI FWD CONTAINER 
COLDPLATE FREON BOY - 4 
COLDPLATE FREON BAY - 5 
COLDPLATE FREON BAY = 6 
COLHPLOTE FREON OUTSIDE FREON BAYS 
TBD 

COLDPLATE FREON MIDSECTION 

FREON PUMP 

TBD 

CABIN AIR. COOLED*.' DIRECT TO HEAT EXCHANGER > 
IMU 

IMU FAN 

NOT ACTIVELY COOLED 

PAYLOAD HEAT EXCHANGER 

COLDPLATE WATER BOY =1 

COLDPLATE WATER BAY =2 

COLDPLATE WATER BAY =30 

COLDPLATE WATER BAY =3B 

CABIN COLDPLATE WATER 

WATER PUMP 

INVERTER =1 

INVERTER =2 

INVERTER =3 

INVER FFR =4 

INVERTFR =5 

FUEL CELL =1 

FUEL CELL =2 

FUEL CELL =3 

FOOD PREP 

HARDCOPY. ENTER ANY CHARACTER TO CONTINUE 


OK1GINAL PAGE 18 
0F POOR QUALITY 


Figure 3.2. 1.2.1 Electrical Power Assignment Menu 
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3*?.J.3 Routine CONSUM 

* * 

The Integration of expended or generated media for the various sources 

assignments Is performed In Routine CONSUM. The control Is through the COMMON/ 

% 

CONS/ which contains the dimensioned variables C0NRAT(20), C0NT0T(20), and 
IAMC0N(20), 


where 


IAMCON(I) » Type of media for source I Indexed as shown on 
Table 3. 2. 1.3.1. 

CONRAT(I) « Media expended (+) or generated (-) for source I 
during most recent time step. 


CONTOT(I) » Total (Integrated) expended (+) or generated for 
source I at current time. 

The source assignment I and the type of media is assigned through the cross 

* 

reference data for the argumented component performance routlne(s) during active 
execution. The arrays CONRAT and CONTOT reflect the algebraic sum of all trans- 
action referencing that source. 

In the latter stages of active execution the IAMCON array Is searched for 
nonzero values. The existence of one or more such values brings up the Initial 
loading Interactive display shown on Figure 3.2. 1.3.1. The COMMON/CMORE/ which 
contains the variable NCON, and the dimensioned variables IC0N(20) and C0NSTA(20) 
Is loaded at this time, 


where 


NCON * Number of nonzero elements In IAMCON, 

ICON » Contains the NCON referenced source numbers In the 
first NCON locations, 

3.2. 1.3.1 


6 

and 

CONSTA(I) ■ Initial quantity for source I. 

It should be noted that the ICON array Is used for Initial loading and 
output (see Routine CONPR) control rather than a repeated search of IANGON. The 
ICON array Includes only those sources referenced Immediately prior to Initial 
loading. Accordingly » any dynamic communication must refer only to sources 
that have been referenced at this time. Additional source references may be 
enforced by setting 

IAMCON(I) > J. 

where I Is the source number and J Is the media Index. Such an entry Is best 
affected Immediately after the call to START and outside the timing loop. 

0 Quantity remaining for each of the sources Is calculated In the output 

Routine CONPR. 
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0 


Table 3. 2. 1.3.1. Consumables Media Index for Source I 


IAMCON(I) 


SYMBOL 

UNITS 

1 

Potable water 

h 2 o 

LSB 

2 

Water 

h 2 o 

LBS 

3 

Carbon Dioxide 

C0 2 

LBS 

4 

Oxygen 

°2 

LBS 

5 

Hydrogen 

H 2 

LBS 

6 

Nitrogen 

n 2 

LBS 

7 

Methane 

ch 4 

LBS 

8 

Hydrogen Peroxide 

H 2°2 

LBS 

9 

Ammonia 

nh 3 

LBS 

10 

Other 

m 

LBS 

11 

Electrical Power 

- 

WATT HI 

12 

Lithium Hydroxide 

LIOH 

LBS 



\ 

! 

i 


3. 2.1. 3.3 
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ITEM 

SOURCE 

NO. 

CONSUMABLES 
SOURCES 
TYPE OF, 
CONSUMABLE 

UNIT 

INITIAL 

LOAOING 

1 

1 

NITROGEN 

LBS 

.000 

2 

2 

OXYGEN 

LBS 

.000 

3 

3 

LITHIUM HYDROXIDE 

LBS 

.000 

H 

H 

MATER 

LBS 

.000 

6 

6 

ELECTRIC POMER 

MATT HRS 

.000 

6 

6 

POTABLE MATER 

LBS 

.000 


♦ *•*♦♦♦!>*•*♦❖♦♦♦ i|ti|e9jc ♦ ^ # *I«*|«i|*jKi^ ♦♦♦♦at* i®M|n|8i|t>|ti|« «•!♦♦♦# ♦♦<*♦♦<**•*♦♦♦♦>•*♦ ♦♦♦♦♦♦<*•►♦ 
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Figure 3. 2. 1.3.1. Typical Interactive Source Initial 

Loading Quantity Display 


3.2. 1.3.4 


3.2. 1.4 Routine TRAJ 


Reference to component simulation routines which Imply Incident (orbital) 
heating as a boundary condition automatically brings Routine TRAJ Into execution. 
There ere three components of Incident heat resulting from orbital operation. 

(1) Solar radiation directly from the Sun, 

(2) Thermal radiation from the planet being orbited, And 

(3) Albedo (reflected solar radiation) from the planet. 

Routine TRAJ controls the calculations for these three heating components. 

The items required to evaluate the incident heating are: 

(1) Location cf the panel with respect to the vehicle coordinate system, 

(2) Attitude of the spacecraft with respect to a geocentric Inertial or 
local vertical, 

(3) Position of the Sun with respect to the geocentric system, and 

(4) Position of the spacecraft with respect to the geocentric system. 

The first Item Is characteristic of the subject panel. These parameters 

are entered through the subject component simulation routine. Two options are 
available to establish the latter Items. The data may be read In from a 
previously generated trajectory tape or calculated with respect to a prescribed 
set of orbital parameters. Both options use the following coordinate systems. 

(1) The Geocentrlr Inertial System (GCI). The coordinate system origin 
Is at the center of the Earth. The X-axis lies In the equatorial 
plane and points toward the vernal equinox. The Z-axIs passes through 
the North Pole. The Y-axis lies In the equatorial plane and forms a 
right-handed system. 

(2) The Vehicle System (VS). This Is the principal system In which the 
flat plate vehicle geometry Is defined. It differs from the geocentric 

Inertial (GCI) system by the amount of pitch, yaw, and roll. 

3.2. 1.4.1 


The Interactive comrnml cation Is Initially through a console display 
requesting the selection of these options as shown on Figure 3.2. 1.4.1. Item 2 
on this display Is a control parameter related to the dimensional units on the 
tape and Is active only for the tape read option. 

The format of the trajectory tape Is shown on Table 3. 2. 1.4.1. With this 
option In effect dynamic communication should be limited to the characteristic 
of the subject panel. 

Interactive comnunl cation Is further extended through the console display 
shown on Figure 3. 2. 1.4.2 for the calculated trajectory option. The orbited 
parameters are Initially loaded with a default Inertial hold circular equatorial 
orbit with the Sun located out the X-axis. 

The attitude hold key establishes whether the spacecraft defined pitch* 
yaw, and roll are referenced to the inertial (GCI) or the orbital plane (Local 
Vertical System). In the Local Vertical System the X-axis Is along the planet 
to spacecraft vector, the Z-axIs Is perpendicular to the orbital plane, and the 
Y-axis completes a right-handed system. 

The Euler angles which are defined as: 

4 * Rotation about the Z-axIs 

e Rotation about the Y-axis 

t Rotation about the Z-axIs 

and are Illustrated on Figure 3. 2. 1.4. 3. 

The orbital parameters, which are used to calculate the time dependent 
coordinate location of the spacecraft are Illustrated on Table 3. 2. 4. 1.2. 

Dynamic cross reference to these parameters Is shown on Table 3. 2. 4. 1.3. 

Absorbed heat Is processed and assembled for the panels through a call to 
Routine SHAD directly from the component performance routine. (See Routine SHAD. 


3.2. 1.4.2 


Having established the basic parameters the vector manipulation Is perform?* 
to obtain the required angles and the components of Incident heating are 
evaluated as follows: 

Solar Radiation to a Flat Plat e 

The direct solar radiation to a flat plate Is given by 
q s • S cos e 3 


where 

q s ■ Incident energy from the Sun (BTU/HR FT 2 ) 

S ■ Solar constant (BTU/HR FT 2 ) 

e 1 ■ The angle between the vehl cl e-Sun vector and a normal to a 
flat plate element (DEG) 

Planetary Thermal Radiation to a Flat Plat e 

The planetary therval radiation from an Isothermal body to a flat plate Is 
given by: 

«t ■ F t 't 

where 

2 

q t * Incident planetary thermal radiation (BTU/HR FT ) 

F t ■ View factor for a flat plate (nd) 

* Total energy emitted from planet unit area (BTU/HR FT ) 

and 

A „ 

! t • (1 - a) S (1 - a) S/4 


3.2. 1.4.3 









C 


* 



whort S 1$ the solar constant, a Is the planetary albedo constant, Ap Is the 
projected area of the planet, and A^ Is the total surface area of the planet. 
For a sphere, Ap//^ ■ 1/4. 

The view factor F t for a flat plate Is defined by the following: 



- 2 sin 



1 Y.t/IJ - h ,,n [ 2 » 1 "’ 1 ♦ 

ifu * in [ 2 s,n '' *2 


or If 


®2 < 


tan' 


-1 


R 


the expression for F t Is 


f 


cos e« 


where 


H ■ R ♦ h 

R * The radius of the planet (n ml) 
h * Vehicle altitude (n ml) 

e. * The angle between a normal to the flat plate and the radius 
c vector from the planet to the vehicle 


Figure 3. 2. 1.4. 4 shows the geometric relationships for planetary thermal 
radiation from an Isothermal body to a flat plate. 


3.2. 1.4.4 


Plant tarv Albedo to a Flat Plata 


Tim planetary albedo radiation to a flat plate Is given by 
utters 

q # ■ Incident planetary albedo (8TU/KR FT*) 

F # * View factor for a flat plate (ml) 

S ■ Solar constant (BTU/HP FT 2 ) 
a ■ Albedo constant 

The view factor F fi for a flat plate was obtained from Reference 2. This expres- 
sion Is a modified solution for planetary thermal radiation multiplied with a 
correction term, which, since planetary albedo obeys Lambert's Law, accounts for 
the cosine distribution not only with respect to the angular radiation from a 
given area but over the sunlit surface of the planet as well. 

Urn view factor for planetary albedo Is expressed as 

F # • F t (e 2 ,h) ^0.86 * 0.14e“ 0,757 h/R ) 

cos Jo, - |o.l369 (* - e 2 ) 3 cos (e c ) (l - e" 5 * 66 h/R * e 2 )2 )]| 

where 

h • Vehicle altitude (n ml) 

e 2 • Angle between a normal to the flat plate and the radius vector 

e. > The angle between the planet-Sun vector and the radius vector 
1 from the planet to the vehicle 
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y 1. 


0 _ « The angle of rotation of a normal to a flat plate element, , 
c measured from a plane containing ter. planet-Sun vector and the 
radius vector from the planet to the vehicle 


Figure 3.2. 1.4.5 shows the geometric relationships between the parameters. 

When the spacecraft is shadowed by tee planet it is orbiting, the solar 
heating rate to each flat element is set to zero. The position of the space- 
craft relative to the Sun and tee planet Is checked in this subroutine to see 
if occulation has occurred. 
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ORBITAL HEATING 
CONTROL PARAMETERS 

ITEM VALUE UNIT 

*1 CONTROL INDICATOR 2 INTEGER 

1 = READ TAPE 

2 s CALCULATE TRAJECTORY 

* 2 UNIT CONVERSION FOR TAPE 1 INTEGER 

1 = EARTH RADI HER. > 

2 = KILOMETERS (KM. > 

I* ^ * $ * * ❖ ^ ❖ * * ** ^ * Jfc ^ ❖;* ^ * * * *’ ^ * * * * ^ $ $ ❖ * ** * ^ ifc * $ $ * Sfc Jfc** * * * ❖** * 


Figure 3. 2. 1.4.1. Typical Interactive Orbital Heating Control Display 
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Table 3.2.1 .4.1. Edited Trajectory Tape Format 


e 



FORTRAN Tape Number 13 


Identification: Edited Trajectory Information Tape 

Type: Binary 

Density: 800 

1108 File Letter: K 




Word Number Type 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 
19 


R 

R 

R 

R 

R 

R 

R 


20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 


R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 


Record Description 


Description 


Ground Elapsed Time (GET) in decimal hours 

Year 

Month 


Calendar date 


Day 
Hour 
Minute 
Second 

Revolution Counter 

Physical Eclipse Indicator (0 Sun, 1 shadow) 
Noon Indicator (1 noon, 0 other) 

Not Used 
Not Used 

Vehicle vector X 

Vehicle vector Y 

Vehicle vector Z 

Vehicle vector X (x dot) 

Vehicle vector Y (y dot) 

Vehicle vector Z (z dot) 

Vehicle vector R (radius vector to 

spacecraft) 


Direction Cosines XTX 
Direction Cosines XTY 
Direction Cosines XTZ 
Direction Cosines YTX 
Direction Cosines YTY 
Direction Cosines YTZ 
Direction Cosines ZTX 
Direction Cosines ZTY 
Direction Cosines ZTZ 
Sun vectors X SUN 

Sun vectors Y SUN 

Sun vectors Z SUN 

Solar Incidence Angle Beta 
Orbital period 


Units 

HR 


ER or KM 
ER or KM 
ER or KM 
ER/HR or KM/HR 
ER/HR or KM/HR 
ER/HR or KM/HR 

ER or KM 


ER or KM 
ER or KM 
ER or KM 
DEG 
HR 


3.2. 1.4. 8 
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48* ❖ H* *» * * *J* * ifc % * * #t #$ *j* # $ # 4* ■$ <* ♦ <« $ & # * * # *> # H* i$i*$* t|t * if* # ♦ $ * t* si* * *9* *S *1* ** # * # # 


ORBITAL PARAMETERS 


ITEM 


VALUE 

UNIT 

* l 

COMP FREQUENCY 

1 

INTEGER 

• 

2 

ATTITUDE HOLD KEY 

1 s INERTIAL 

2 = LOCAL VERTICAL 

* 

% 

INTEGER 

t. * 

'3 

H 

6 

SUN COORDINATE X . > 

Y 

Z 

2. 346*000* OH * 

o. eeooooo 
0.0000000 

ER. 

ER. 

ER. 

* 

6 

? 

8 

EULER ANGLE ABOUT Z 

Y 

X 

, .000 
.000 
.000 

RAD. 

RAD. 

RAD. 

*9 

OR6IT SEMIMAJOR AXIS 

1.029 

ER. 

io 

ORBIT ECCENTRICITY 

. ooo 

14/D 

ii 

ORBIT INCLINATION 

.000 

RAD. 

12 

RIGHT ASCENSION 

.000 

RAD. 

13 

ARGUMENT OF PERIGEE 

.000 

RAD. 

14 

TIME OF PERIGEE PASSAGE 

.000 

ER. 


* *$c * * * * Jj* * ^ >!< *: * *$< ** * *Ji * H< * H« t* * »*<* * *f * #* •!< * *{« ^ * ** # H» >•* ft ft H« * H< K* * *fe H* * »}<#$« * »6« *8< »|»a®c * # 


Figure 3. 2. 1.4. 2. Typical Interactive Orbital 

Parameter Display 
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Hi" 

V SC axis Initially aligned with OCX axis 



2. Rotate 90 degrees about Z-axIs (Positive rotation Is X-axis towards Y-axis) 



Z 



v ■ 90 DEG 
e • 0 DEG 
0*0 0E6 


3. Rotate 90 degrees about Y-axis (Positive rotation Is Z-ax1s towards X-axis) 


4. 



Rotate 90 degrees about X-axis (Positive rotation is Y-axis towards Z-axIs) 




V • 90 DEG 
0 ■ -90 DEG 
0 • -90 DEG 


Figure 3. 2. 1.4. 3. Illustration of Euler Angle Sequence 

Showing Spacecraft Position 


3.2.1.4.10 
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Table 3.2.1 .4.2. Classical Elements of an Orbit 



c 

The elements which define an orbit are: 
a - Semi major axis of the orbit 
e - Eccentricity of the orbit 

1 - Inclination of the orbit, angle between the orbital 
plane and the equator 0° <_ 1 <, 180° 

ft - Right ascension of the ascending node 0° <: si < 360° 

a) - Argument of perigee 

x - Time of perigee passage 

i 

! 

i 

i 

i 

G 

i 
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3.2.1. 5 Routine SHAD 


The absorbed heat for a panel 1 Is assembled In Routine SHAD In response 
to a direct call from the component routine processing node 1. 

If the component routine processing node 1 does not reference shadowing the 
assembly Is simply 

<»ABS " ° (< *s + < *a ) * e q t 

from the three Incident heat components calculated by control throuph Routine TRAJ 
and the panel solar absorbtlvlty and thermal emlsslvlty prescribed for node 1. 

If the component routine processing node 1 references a shadowing node j 
and a stand-off vector storage location n the shadowing effect Is taken into 
account before the absorbed heat Is calculated. 

The shadowing calculations set up two concentric solid angles defined by 
the angles of lune and a 2 . which define the total and partial shadowing 
bounds for solar radiation as shown for a simple configuration on Figure 3. 2. 1.5.1. 
For simplicity the configuration Is representative of a system In which the 
normal to 1 and j and the stand-off vector coincide. The panels mgy have any 
relative position on the spacecraft defined by their dihedral angle (b) and 
angle of Incidence (a). Slmularly for the stand-off vector parameters B m and 
o^. The parameters u m and B m are the angle of Incidence and dihedral angles 
of a panel that result In a normal that Is In the direction of a vector drawn 
from the center of panel 1 to the center of panel j. Panels 1 and j are 
separated by the distance m along this vector. Albedo and thermal shadowing 
effects are evaluated as: 

°<q, -q,, -q^ 


3.2. 1.5.1 



and 

m *ti * q tj 

for node 1. 

Interactive comnunl cation with the shadowing node and stand-off vector 
parameters Is through the console display shown on Figure 3.1. 2.6.3. 

The cross reference to shadowed node, shadowing node and stand-off vector 
parameters are shown on Table 3.2.1. 5.1. 
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3.2. 1.5.2 


AllfiLE tfETNEEN SUN VECTOR AND STAND OFF VECTOR 



Figure 3.2 1.5.1. Typical Solar Shadowing Illustration 


3.2. 1.5.3 






Table 3.2. 1.5.1. Dynamic Cross Reference Parameters 
for Node I Shadowed by Node j With 
Stand-Off Vector Storage In N 


VARIABLE 

SYMBOL 

Shadowed node area 

A t 

Shadowing node area 

A J 

Angle of Incidence shadowed 
node 

a 1 

Angle of Incidence shadowing 
node 

a J 

Angle of incidence stand-off 
vector (equivalent) 

a m 

Dihedral angle shadowed node 


Dihedral angle shadowing node 


Dihedral angle stand-off 
vector (equivalent) 

e m 

Stand-off distance 

m 


TTTi] 


C(J,10) 

C(2,11) 

C(J,11) 

c(m,ii) 

C(1.12) 

C(J,12) 

C(M,12) 

C(M,10) 
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3.2*2 Output Routines 

The basic control of output date Is through the control Routine PRINT via 
Initial Interactive Print Option Display shown on Figure 3. 1.1. 2.2. This 
s^tlon discusses the subsequent Interactive, and possible dynamic, communl cation 
resulting from the selection ot various print options, os well as various output 
functions which are automatically executed In response to the particular model 
component performance routines characteristics. A summary of the output 
routines Is shown on Table 3.2.2. 1. 


3.2.2. l 



Table 3. 2.2.1. Summary of Output Routines 

\ • • 


NODPKT Displays fluid properly data with generic names 

at selected print frequency. 

NASPRT Displays fluid property data with assigned node 

names at selected print frequency. * 


6ASPR Displays atmospheric property data at selected 

print frequency. 

CONPRT Displays consumables data for referenced sources 

at selected print frequency. 

PLOOT Stores data for plotting. 

SCHEM Displays modelled system schematic. 
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O 3.2.2.I. Routine NOOPRT 

The basic output of fluid property data Is written from Routine NOOPRT. 
The routine automatically outputs data during passive execution for all 
referenced node numbers in the order of reference as shown on Figure 3.2.2.1.1 
unless the user specifies the select node option* The output at the print 
frequency specified at the time of Initialization control (Figure 3.1. 1.1*3). 

u I ^ 

The generic name of the referenced, nodes Is Included. 


O 


ft 


3.2.2. 1.1 
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TFEAR 

TEST TFEAR PLOT OPTION 


FLUID PROPERTIES 


TIME 

NODE. 

= 1.500 


COMP 

FLUID 

WCP 

HEAT 



TEMP 

TEMP 

BTU/HR 

LOAD 

NO 



DEG 

DEG 

DEG 

BTU/HR 

l 

CA3IN IN 

✓HEATER OUT 


498.364 

501.000 

2001.000 

2 

CABIN OUT 

✓LIOH IN 


503.450 

501.000 


3 

LI OH OUT 

✓COND. IN 


503.870 

501.000 


4 

COND. OUT 

/HEATER IN 

498. 448 

496.573 

501.000 

450.000 

5 

EVAP IN 

✓JUNCTION 

495. OOO 

590.630 

1490.483 


6 

EVAP OUT 

✓PLATE IN 

495.601 

495.755 

1490.483 


? 

PLATE OUT 

✓INTF IN 


495. 625 

1490.483 

3793. 755 

8 

INTF OUT 

✓PLATE IN 

499. 816 

498. 170 

1490.483 

150. 000 

9 ■ 

PLATE OUT 

✓EXCH IN 


499.598 

1490.403 


10 

EXCH OUT 

✓DIVERT 


501.023 

1490.483 


! 12 

LEG 

✓BRANCH 


501.023 

7.452 


11 

LEG 

✓MIXER 


501.023 

1483.031 


! 13 

LEG 

✓RADIN 

482.899 

501.023 

3.726 


1 17 

LEG 

✓RADIN 

482. 899 

501.023 

3.726 


14 

RADOUT 

✓RADIN 

475. 798 

483. 158 

3.726 


! ( ) IS 

RADOUT 

✓RADIN 

470.936 

476.049 

3.726 


1 V - 16 

RADOUT 

✓MIXER 


471 . 178 

3.726 


! 18 

RADOUT 

✓RADIN 

475. 798 

483. 158 

3. 726 


I 19 

RADOUT 

✓RADIN 

371 . 416 

476.049 

3.726 


i 20 

RADOUT 

✓MIXER 


373 . 764 

3.?26 


i 21 

EXCH IN 

✓PLATE OUT 


522.508 

100.000 


1 2? 

EXCH OUT 

✓PLATE IN 

523.615 

501 . 291 

100. 000 

200.000 
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Figure 3 . 2 . 2 . 1 . 1 . Typical Fluid Property Output Data - Generic Names 
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3. 2*2. 2 Routine NASPRT 


If the option to select nodes In Is effect* output Interactive communication 
Is extended through a display requesting the node numbers to be output and 
whether or not to use assigned or generic node names. If the generic name option 
Is In effect* the output ts through Routine NODPRT, except that the selected set 
of node numbers (In the order they were Input)* rather than all referenced model 
nodes are printed. Selection of the option of assigning node names extends the 
Interactive communication through a display requesting the nodal names In tee 
same order as the node number selection. Subsequent output Is through Routine 
NASPRT reflecting the selected node numbers and the assigned names as shown on 
Figure 3. 2. 2. 2.1. 
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TFEfiR 

TEST TOPE REOQ CAPABILITY 

FLUID PROPERTIES 

.000 

COMP 
TEMP 
DEG 

I'flK r hi 

FUEL CELL 521. OOO 

EQUIPMENT 521.000 

POYL.OOD 521.000 


HEOT 

LORD 

BTU/HR 


TIME = 

NOOE 
NO 

«L 
6 
8 

22 
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FLUID 

TEMP 

DEG 

510. 165 
498.474 
515.559 
S2Q.339 


WCP 

BTU/HR 

DEG 

5O1.OO0 
1490.484 
1490.484 
100. OPO 


Figure 3. 2. 2. 2.1. Typical Fluid Property Output 

Data - Assigned Names 
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3.2.2. 3 Routine GASPR 


Reference to atmospheric properties for one or more component performance 
routines will automatically execute the atmospheric data output Routine GASPR. 
The data as shown on Figure 3.2. 2.3.1 Is at the selected print frequency and 
Includes all node numbers referencing atmospheric processing. 




GRS PROPERTIES 


TIME = 1 . 50Q 


NODE 

SPECIFIC 

WATER 

NITROGEN 

OXYGEN 

C02 

TOTRL 

NO 

HERT 

PSI 

PSI 

PSI 

PSI 

PSI 

1 

.210 

. 106 

11.600 

3.100 

.073 

14.700 

2 

.210 

. 10? 

11.600 

3. 100 

.073 

14. 700 

3 

.210 

. 107 

11.600 

3. 100 

.073 

14.700 

4 

.210 

. 106 

11.600 

3. 100 

.073 

14.700 
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Figure 3. 2. 2. 3.1. Typical Atmospheric Data Output 
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3. 2. 2. 4 Routine CONPRT 


Reference to one or more storage sources will automatically execute Routine 
CONPRT to output the summary of the status of the various sources. The output 
Is at the selected print frequency and reflects the remaining quantity at the 
given time as shown on Figure 3. 2. 2.4.1. The remaining quantity for a source 1 
Is calculated Internally as 

CONSTA(I) - CONTOT(I) 

referenced In the discussion of Routine CONSUM. Dynamic communication with the 
remaining quantity should reference these variables. 
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CONSUMABLES 

USAGE 


TIME = 

1.500 

INITIAL 

QUANTITY 

QUANTITY 

SOURCE 

AVAILABLE 

USED 

REMAINING 

1 

.000 

2.298 

-2.298 

2 

.000 

.848 

-.848 

3 

.000 

.000 

.000 

4 

.000 

-.324 

.324 

S 

.000 

197.947 

-197.947 

6 

.000 

11.900 

-11.900 


C 


ORIGINAL PAGE E 
OP POOR QUALITY 



Figure 3. 2. 2. 4.1. Typical Consumables Oata Output 
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3. 2.2. 5 Routine PLOOT 


Selection of the Plot Option will extend the Interactive communication 
through a display requesting the node numbers and type of Information to be 
plotted. A completed display of this type Is shown on Figure 3.2.2.5. 

The Plot Option automatically executes Routine PLOOT which stores the data 
to be plotted at each calculated time point. Near run termination the data 
stored Is prepared In screen plots by Routine DRAM. 


************************************ 
FUNCTION NENU 
DESCRIPTION 


VPE 

1 

2 

3 

M 

6 

6 

7 

8 
9 

1© 

11 


FlUIO TEMP 
COMP TENP 
FLON RATE 
HEAT 

ACC CONSUN 
HEATER PNR 
PRT PRES H20 
PRT PRES N2 
PRT PRES 02 
PRT PRES C02 
HTR ENERGY 


( DEG/R > 
< DEG/R > 
(LBS/HR) 
< BTU/HR) 
(LB8> 
(BTU/HR) 
(PSD 
(PSD 
(PSD 
<MMHG) 
(BTU> 


ate * i|t * ifc ift a|i «|t ife $ $ * $ $ $ its H* i|s ifc ifc ifc ifc* ifc $ * ifc * s|e * jft ifc ifc* 


♦ * * * * * * * ** He * * aft * * 4c * * * He * * * * $ * * * * ** * * * * * >|c * * * * * * * * ** * * ate * % * ifc * * 

PLOT CONTROL 

TEN NODE NO TYPE HAX HIN 

1 2 FLUID TENP .000 .088 

2 12 FLOH RATE .000 .008 

3 22 HEAT . 000 . 088 

*********** ********111**41**111*111*****111111****111****************** 


Figure 3.2.2.S. Typical Plot Control Data Display 


3.2. 2. 6 Routine SCHEM 


Request fo* a schematic of the modelled system will automatically execute 
Routine SCHEM Immediately prior to passive execution. Routine SCHEM processes 
the referenced node numbers and their penerlc names to asslpn a node number and 
type to each subplot In a composite display of the system schematic. Routine 
SOCH prepares the control parameters for a single page of the schematics and 
then processes a call to Routine PICT to display that page. The paging 
process continues until the schematic Is complete. A typical single page 
schematic Is shown on Figure 3.2. 2. 6.1. 


3. 2.2. 6.1 
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